On the role of entanglement and correlations in mixed-state quantum computation * 
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In a quantum computation with pure states, the generation of large amounts of entanglement is 
known to be necessary for a speed-up with respect to classical computations. However, examples 
of quantum computations with mixed states are known, such as the DQC1 model [E. Knill and 
R. Laflamme, Phys. Rev. Lett. 81, 5672 (1998)], in which entanglement is at most marginally 
present, and yet a computational speed-up is believed to occur. Correlations, and not entanglement, 
have been identified as a necessary ingredient for mixed-state quantum computation speed-ups. Here 
we show that correlations, as measured through the operator Schmidt rank, are indeed present in 
large amounts in the DQC1 circuit. This provides evidence for the preclusion of efficient classical 
simulation of DQC1 by means of a whole class of classical simulation algorithms, thereby reinforcing 
the conjecture that DQC1 leads to a genuine quantum computational speed-up. 



PACS numbers: 3.67.Lx 



I. INTRODUCTION 



Quantum computation owes its popularity to the real- 
ization, more than a decade ago, that the factorization 
of large numbers can be solved exponentially faster by 
evolving quantum systems than via any known classical 
algorithm [l[. Since then, progress in our understand- 
ing of what makes quantum evolutions computationally 
more powerful than a classical computer has been scarce. 
A step forward, however, was achieved by identifying en- 
tanglement as a necessary resource for quantum compu- 
tational speed-ups. Indeed, a speed-up is only possible if 
in a quantum computation, entanglement spreads over an 
adequately large number of qubits Q. In addition, the 
amount of entanglement, as measured by the Schmidt 
rank of a certain set of bipartitions of the system, needs 
to grow sufficiently with the size of the computation 
Whenever either of these two conditions is not met, the 
quantum evolution can be efficiently simulated on a clas- 
sical computer. These conditions (which are particular 
examples of subsequent, stronger classical simulation re- 
sults based on tree tensor networks (TTN) Q) are only 
necessary, and thus not sufficient, so that the presence 
of large amounts of entanglement spreading over many 
qubits does not guarantee a computational speed-up, as 
exemplified by the Gottesman-Knill theorem [5j|. 

The above results refer exclusively to quantum com- 
putations with pure states. The scenario for mixed-state 
quantum computation is rather different. The intriguing 
deterministic quantum computation with one quantum bit 
(DQC1 or 'the power of one qubit') @ involves a highly 
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mixed state that does not contain much entanglement 
and yet it performs a task, the computation with fixed 
accuracy of the normalized trace of a unitary matrix, 
exponentially faster than any known classical algorithm. 
This also provides an exponential speedup over the best 
known classical algorithm for simulations of some quan- 
tum processes |8[. Thus, in the case of a mixed-state 
quantum computation, a large amount of entanglement 
does not seem to be necessary to obtain a speed-up with 
respect to classical computers. 

A simple, unified explanation for the pure-state and 
mixed-state scenarios is possible Q by noticing that the 
decisive ingredient in both cases is the presence of correla- 
tions. Indeed, let us consider the Schmidt decomposition 
of a vector |\&), given by 
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where {ia\3a) — == $ij and X is the rank of the 

reduced density matrices pa = Tr B[\ip)(ip\] an d Pb = 
Tcj\.[\ip){ip\]j and the (operator) Schmidt decomposition 
of a density matrix p given by [9( 
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parts A and B, with x t = x ' 2 it p = Let the 

density matrix p t denote the evolving state of the quan- 
tum computer during a computation. Notice that p t can 
represent both pure and mixed states. Then, as shown 
in Refs. Q and Q, the quantum computation can be ef- 
ficiently simulated on a classical computer using a TTN 
decomposition if the Schmidt rank x of p according to 
a certain set of bipartitions A : B of the qubits scales 
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polynomially with the size of the computation. In other 
words, a necessary condition for a computational speed- 
up is that correlations, as measured by the Schmidt rank 
x", grow super-polynomially in the number of qubits. In 
the case of pure states (where x — v?) these correla- 
tions are entirely due to entanglement, while for mixed 
states they may be quantum or classical. 

Our endeavor in this paper is to study the DQC1 
model of quantum computation following the above line 
of thought. In particular, we elucidate whether DQC1 
can be efficiently simulated with any classical algorithm, 
such as those in [3, 0] (and, implicitly, in (2|), that ex- 
ploits limits on the amount of correlations, in the sense of 
a small x according to certain bipartitions of the qubits. 
We will argue here that the state pt of a quantum com- 
puter implementing the DQC1 model displays an expo- 
nentially large x i m spite of it containing only a small 
amount of entanglement 0. We will conclude, there- 
fore, that none of the simulation techniques mentioned 
above can be used to efficiently simulate 'the power of 
one qubit'. 

On the one hand, our result indicates that a large 
amount of classical correlations are behind the (sus- 
pected) computational speed-up of DQC1. On the other 
hand, by showing the failure of a whole class of clas- 
sical algorithms to efficiently simulate this mixed-state 
quantum computation, we reinforce the conjecture that 
DQC1 leads indeed to an exponential speed-up. We note, 
however, that our result does not rule out the possibility 
that this circuit could be simulated efficiently using some 
other classical algorithm. 

II. DQC1 AND TREE TENSOR NETWORKS 
(TTN) 



[H 0, 13 require that p t be efficiently represented with a 
TTN [J] (or a more restrictive structure, such as a prod- 
uct of A:-qubit states for fixed k [2:] or a matrix product 
state 0) at all times t G [0, T\. Here we will show that 
the final state pt, henceforth denoted simply by p, can- 
not be efficiently represented with a TTN. This already 
implies that none of the algorithms in 0, 0, 0| can be 
used to efficiently simulate the DQC1 model. 

Storing and manipulating a TTN requires computa- 
tional space and time that grows linearly in the number 
of qubits n and as a small power of its rank q. The rank 
q of a TTN is the maximum Schmidt rank x\ over all 
bipartitions Aj : Bi of the qubits according to a given 
tree graph whose leaves are the qubits of our system. 
See [4( for details. The key observation of this paper 
is that for a typical unitary matrix U n , the density ma- 
trix p in Eq. (|2.2|) is such that any TTN decomposi- 
tion has exponentially large rank q. By typical, here we 
mean a unitary matrix U n efficiently generated through 
a (random) quantum circuit. That is, U n is the product 
of poly(n) one-qubit and two-qubit gates. In the next 
section we present numerical results that unambiguously 
suggest that, indeed, typical U n necessarily lead to TTN 
with exponentially large rank q. 

We notice that the results of the next section do not 
exclude the possibility that the quantum computation in 
the DQC1 model can be efficiently simulated with a TTN 
for particular choices of U„ . For instance, if U n factorizes 
into single-qubit gates, then p can be seen to be efficiently 
represented with a TTN of rank 3, and we can not rule 
out an efficient simulation of the power of one qubit for 
that case. Of course, this is to be expected, given that 
the trace of such U n can be computed efficiently in the 
first place. 



The DQC1 model, represented in Eq. I|2.1[) . provides 
an estimate of the normalized trace Ti(U n )/2 n of a n- 
qubit unitary matrix U n G U(2 n ) with fixed accuracy 
efficiently 6] . For discussions on the classical complexity 
of evaluating the normalized trace of a unitary matrix, 
see 0. 
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This quantum circuit transforms the highly-mixed initial 
state po = |0)(0| ®I n /2 n at time t = into the final state 
Pt at time t = T, 
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through a series of intermediate states pt, t G [0, T\. The 
simulation algorithms relevant in the present discussion 



III. EXPONENTIAL GROWTH OF SCHMIDT 
RANKS 

In this section we study the rank q of any TTN for 
the final state p of the DQC1 circuit, Eq. (|2~2]l . We 
numerically determine that a lower bound to such a rank 
grows exponentially with the number of qubits n. 

The Schmidt rank x of a pure state \p,t> A ip B ) 

\ PMb ) = p \4> a )\^b) = °iA\<f>A)®O iB \i; B ) (3.1) 

i=l 

obtained by applying the density matrix p onto a product 
state \4>a)\4 } b) is a lower bound on the operator Schmidt 
rank x of p, i.e., x" > X- For the purpose of our numer- 
ics, we consider the pure state J7„|0)® n . We build U n as a 
sequence of 2n random two-qubit gates, applied to pairs 
of qubits, also chosen at random. The random two-qubit 
unitaries are generated using the mixing algorithm pre- 
sented in [l(|. Note that applying 2n gates means that 
the resulting unitary is efficiently implementable, a sit- 
uation for which the DQC1 model is valid. For an even 
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number of qubits n, we calculate the smallest Schmidt 
rank \ over all n/2 : n/2 partitions of the qubits (sim- 
ilar results can be obtained for odd n). The resulting 
numbers are plotted in Fig JTJ). 
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FIG. 1: (Color Online) :Lower bound for the operator Schmidt 
rank x of the DQC1 state for any equipartition n/2 : n/2, as 
given by the Schmidt rank \ °f the pure state in Eq. (|3.1[1 . 
The dots are for even numbers of qubits, and the fit is the line 
2 n//2 . x m calculated for a pure state obtained by applying 
2n random 2-qubit gates on the state |0)® n . This is evidence 
that for a typical unitary U n , the rank q of any TTN for the 
DQCI state p in Eq. (|2,2[l grows exponentially with n. 



The above numerical results strongly suggest that the 
final state p in the DQCI circuit has exponential Schmidt 
rank for a typical unitary U n . We are not able to pro- 
vide a formal proof of this fact. This is due to a general 
difficulty in describing properties of the set U 9C (2") of 
unitary matrices that can be efficiently realized through 
a quantum computation. Instead, the discussion is much 
simpler for the set U(2") of generic n-qubit unitary ma- 
trices, where it is possible to prove that p cannot be ef- 
ficiently represented with a TTN for a Haar generated 
U n € U(2"), as discussed in the next section. Notice 
that Ref. [11| claims that random (but efficient) quan- 
tum circuits generate random n-qubit gates U n € U 9C (2") 
according to a measure that converges to the Haar mea- 
sure in 0(2"). Combined with the theorem in the next 
section, this would constitute a formal proof of the other- 
wise numerically evident exponential growth of the rank 
q of any TTN for the DQCI final state p. 



IV. A FORMAL PROOF FOR THE 
HAAR-DISTRIBUTED CASE 

Our objective in this section is to analyze the Schmidt 
rank x" of the density matrix p in Eq. (|2.2p for certain 



bipartitions of the n+1 qubits, assuming that U n € U(2") 
is Haar-distributed. 

It is not difficult to deduce that for any tree of the n+1 
qubits, there exists at least one edge that splits the tree 
in two parts A and B, with ua and ng qubits, where 
no = min(n,4,«s) fulfills n/5 < no < 2n/5. In other 
words, if a rank-g TTN exists for the p in Eq. (|2.2p . then 
there is a bipartition of the n + I qubits with no qubits 
on either A or B and such that the Schmidt rank x < <7- 
Theorem [TJ our main technical result, shows that if U n is 
chosen randomly according to the Haar measure, then the 
Schmidt rank of any such bipartition fulfills x" > O(2 no ). 
Therefore for a randomly generated U n £ U(2 n ), a TTN 
for p has rank q (and computational cost) exponential in 
n, and none of the techniques of 0, d, 0] can simulate the 
outcome of the DQCI model efficiently. 

Consider now any bipartition A : B of the n+1 qubits, 
where A and B contain ua and ns qubits, with the min- 
imum no of those restricted by n/5 < no < 2n/5. With- 
out loss of generality we can assume that the top qubit 
lies in A. Actually, we can also assume that A contains 
the top via qubits. Indeed, suppose A does not have the 
via top qubits. Then we can use a permutation P n on all 
the n qubits to bring the tla qubits of A to the top tla 
positions. This will certainly modify p, but since 
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where V n = P n U n Pl is another Haar-distributed uni- 
tary, we obtain that the new density matrix is of the 
same form as p. Finally, in order to ease the notation, 
we will assume that ua — no (identical results can be 
derived for ub — no). Thus n/5 < ua < 2n/5. 
We note that 
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so that if we multiply p by the product state 
\(t>s) = \t,i,j) = 



®U n , 
(4.2) 

(4.3) 



where a = (t,i, j), t — 0, 1; i — 1, . . . <1a] j = 1, • • • ds, 
we obtain \ips) = p\4>s) where 

U \_ / zkr{%h3) + \V®U n \i,3)) if t = 
Wa) \ 3 ^ r (|i,i J j) + |0)(8^|i > i» if * = 1 l4 - 4j 

This also justifies our choice of the pure state used in the 
numerical calculations in the previous section. 
Let us consider now the reduced density matrix 



Tr A [|Vs)(^|] 

(|j><j| + Tr A [U n \i,j)(i,j\Ul}) (4.5) 



for t = (for t — 1, U n and need to be exchanged). 
For a unitary matrix U n randomly chosen according to 
the Haar measure on U(n), U n \i,j) is a random pure 
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state on A® B. Here, and henceforth A is the space of 
the first riA qubits without the top qubit. It follows from 
[lH that the operator 



Q = Tr A [U n \i,j){i,j\Ul] 



(4.6) 



has rank d A ■ Therefore the rank of a? (equivalently, the 
Schmidt rank x of \ip s )) is at least 2™°. From Eq. ([XT]) 
we conclude that the Schmidt rank of p fulfills > 2™° > 
2™/ 5 . We can now collate these results into 

Theorem 1 Let U n be an n-qubit unitary transforma- 
tion chosen randomly according to the Haar measure on 
U(2"), and let A: B denote a bipartition of n + 1 qubits 
into riA and ub qubits, where no = min(?i J 4, ns)- Then 
n/5 < no < 2n/5 and the Schmidt decomposition of p in 
Eq. \2.ty according to bipartition A : B fulfills x S > 2 n/5 . 



We have seen that we cannot efficiently simulate DQC1 
with an algorithm that relies on having a TTN for p with 
low rank q. However, in order to make this result robust, 
we need to also show that p canot be well approximated 
by another p accepting an efficient TTN. We do this in 
Appendix O 



V. CONCLUSIONS 

The results in this paper show that the algorithms of 
are unable to efficiently simulate a DQC1 circuit. 
The efficiency of a quantum simulation using these algo- 
rithms relies on the possibility of efficiently decomposing 
the state p of the quantum computer using a TTN. We 
have seen that for the final state of the DQC1 circuit no 
efficient TTN exists. 

It is also interesting to note that the numerics and 
Theorems 1 and 2 in this paper can be generalized for 
any fixed polarization r, (0 < r < 1) of the initial state 
t|0)(0| + (1 - t)I/2 of the top qubit of the circuit in Eq 
(|2.ip , implying that the algorithms of @j S 13] are also 
unable to efficiently simulate the power of even the tiniest 
fraction of a qubit. 
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APPENDIX A: DISTRIBUTION OF THE 
SCHMIDT COEFFICIENTS 

In this Appendix we explore the robustness of the 
statement of Theorem [JJ To this end, we consider the 



Schmidt rank \ for a density matrix p that approximates 
p according to a fidelity F(0\ , O2) defined in terms of the 
natural inner product on the space of linear operators, 




Tr^OxWTr^Oa), 



F(0 1) 2 ) = Tr(0\0 



where F = 1 if and only if 0\ = O2 and F = \(ipi\ip2)\ 2 
for projectors Oi = on pure states \ipi). We will 
show that if p is close to p, then x for a bipartition as 
in Theorem I is also exponential. To prove this, we will 
require a few lemmas which we now present. 

Lemma 1 Let |\&) be a bipartite vector with x terms in 
its Schmidt decomposition, 

x x 
|*> = JV*5^Ai|iA>|iB>, Xi>K+i>0, X)A? = 1, 



where Nq, = \/(W\W), and let |<I>) be a bipartite vector 
with norm 7V$ and Schmidt rank x' , where x' <X- Then, 



max|(*|$}| = NyN^t 



x 



(Al) 



Proof: Let p,i denote the Schmidt coefficients of |$). 
It follows from Lemma 1 in [l^] that max|$) KvE"^}! = 

N-q,N<£ XiPi, and the maximization over pi is done 
next. A straightforward application of the method of 
Lagrange multipliers provides us with pi = c\i, i = 
1,2, ...,x' f° r some constant c. Since 52iLiM? = 1 = 
^x' 



Tii^, C =l/ v ELiA?.Thus, 



max|(*|$)| = cA^A^^A^ 

i=l 

and the result follows. □ 
We will also use two basic results related to majoriza- 
tion theory. Recall that, by definition, a decreasingly or- 
dered probability distribution p — (jpi,p2, ■ ■ ■ ,Pd), where 
Pa > Pa+i > 0, J2 a Pa = 1; is majorized by another such 
probability distribution q, denoted p -< q, if q is more 
ordered or concentrated than p (equivalently, p is flatter 
or more mixed than q) in the sense that the following 
inequalities are fulfilled: 



^2p a <^2q a Vfc = l,. 



(A2) 



a— 1 a— 1 



with equality for k = d. The following result can be 
found in Exercise II. 1.1 5 of [l^ : 

Lemma 2 Let pg and p$ be density matrices with eigen- 
values given by probability distributions x and y. Let 
o~(M) denote the decreasingly ordered eigenvalues of her - 
mitian operator M . Then 

o-{p3 + Py) <x + y- 
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The next result follows by direct inspection. 



Proof: For any product vector of Eq. (|4. 3|) we have 



Lemma 3 Let coefficients Si, 1 < i < d, be such that 
—8 < Si < 5 for some positive 5 < 1 and 8% = 1, and 
consider the probability distribution p({5i}), 



2d 



2d 



l + S 
2d 



Then 



where 



6* = 



5 i<d/2 
-S i > d/2 



and we assume d to be even. 

Finally, we need a result from [l3T ] : 
Lemma 4 With probability very close to 1, 



Pr 



(1_ 5 )^L<Q< (1 + 5)1- 



> 1 



10 d A 



2d a 



■)(-d B 5 2 /141n2) 



>l-0 



exp(<5 2 cxp(n)) 



(A3) 



where d A = 2 nA = 2"° and d B = 2 nB = 2"-™ n+1 , and the 
operator Q defined in Eq. |^.6'| ) is within a ball of radius 
8 of a (unnormalized) projector T/dA of rank dA [pro- 
vided ds is a large multiple of dA^ogdA/S 2 Jl$i] , which 
is satisfied for large n, given that n/5 < no < 2n/5]. 

Our second theorem uses the fact that the Schmidt de- 
composition of p does not only have exponentially many 
coefficients, but that these are roughly of the same size. 

Theorem 2 Let p, U n , and A:B be defined as in Theo- 
rem^ If F(p,p) > 1 — e, then with probability p(S,n) = 
1 — 0(exp(— S 2 exp(rt))), the Schmidt rank for p according 
to bipartition A:B satisfies x" > (1 — 4e — <5)2™/ 5 . 



\(tij\pp\tij)\ < N S N S 

< N s N s g{^/d A ) 



\ k=l 



(A4) 



where 



1 + (1 + S)x 



(A5) 



and N s = vWFN), ^5 = VWWWjj- The first 
inequality in (|A4[) follows from Lemma 1, whereas the 
second one follows from the fact that the spectrum p of 



PB = (N s )- 2 TT A [p\Uj)(Uj\p} = + Q), 

where Q has all its d A non-zero eigenvalues qi in the 
interval 2 _ ™°(1 - 8) < q t < 2 _ ™°(1 +8), is majorized by 
p({6*}), as follows from Lemmas 2 and 3. Then, 



l-e< 



Tipp 



y/Trp 2 y/Trp 2 

J2s ("Iwl") 



vS («'IP 2 I«')E 5 - (a"|p 2 |«") 



< 9(x*/d A ), 



E S '(N 3 ,) 2 J: S 4N S ») 2 



where in the last step we have used the Cauchy-Schwarz 
inequality, |(x|y)| < y/ (x\x) \J (y\y). The result of the 
theorem follows from g(x^/2 n °) > 1 — e. ^ 
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